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Introduction 

This  final  report  describes  research  completed  under 
Contract  AP49(638)-1243,  during  the  period  1  July  1963  to 
31  December,  1969*  The  primary  objectl/es  of  the  project 
were  as  follows: 

(1)  The  development  of  theoretical  methods  for 
treating  the  amplitude  and  phase  characteristics  of  body 
waves  In  the  earth. 

(2)  The  development  of  parallel  numerical  methods  for 
computing  the  amplitudes  and  phases  of  body  waves  In  a 
radially  heterogeneous  earth, 

(3)  The  application  of  these  results  to  a  more  com¬ 
plete  understanding  of  observed  long-  and  short-period 
signals  from  explosions  and  earthquakes. 

Major  results: 

The  primary  outcome  of  the  research  may  be  summarized* 

(1)  Further  understanding  of  the  effect  of  the  crust 
on  long-period  body  waves.  In  Phlnney  (1964)  and  Phlnney 
(1965)1  the  distortion  effect  of  the  crust  on  teleselsmlc  P 
waves,  and  the  theory  of  the  body  wave  portion  of  the  crustal 
long-period  seismogram  were  treated, 

(2)  A  full  treatment  of  the  effects  produced  by  the 
core  boundary  and  the  lower  mantle  on  body  waves  Is  now 
available,  from  a  series  of  papers  beginning  In  1966,  The 
core  diffraction  proved  to  be  especially  Important  for 
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testing  these  ideas  because  it  is  about  the  highest-quality, 
least>oluttered  signal  on  exlstintr  roc^ris. 

(3)  A  full  theory  for  the  interaction  of  body  waves 
with  regions  of  high  velocity  gradients  and  with  negative 
gradients  is  now  available  (Phlnney,  1970  and  Chapman,  1969). 

The  kinematic  behavior  of  waves  in  a  shadow  region  may  be 
described. 

(4)  The  development  of  a  scheme  for  the  exact  numeri¬ 
cal  representation  of  a  long-period  body  wave  signal.  This 
scheme  has  been  Implemented  in  an  existing  computer  program, 
and  is  being  applied  in  the  synthesis  of  signals  from  a 
realistic  earth  model. 

(5)  The  discovery  of  the  utility  of  a  complex  fre¬ 
quency  variable  as  a  means  of  smoothing  data,  with  maximum 
weight  given  to  early-arriving  signals. 

Most  of  the  published  theoretical  and  numerical 
results,  in  addition  to  a  great  deal  more  in  the  form  of 
notes,  represents  the  source  material  for  a  monograph  now  in 
preparation  by  the  Principal  Investigator.  The  publication 
of  this  material  will  probably  take  another  18  months. 

Resources  generated; 

In  addition  to  the  principal  published  results,  the 
following  resources  for  data  analysis  and  theoretical  work 
were  generated: 

(1)  Computer  programs  for  (a)  calculation  of  the 

response  of  a  layered  crust  to  body  waves;  (b)  cal- 
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culation  of  the  signal  generated  by  a  point  source 
in  a  layered  medium;  (c)  calculation  of  the  signal 
generated  by  a  point  source  in  a  spherical  earth; 

(d)  reduction  and  analysis  of  long-period  core 
diffractions;  (e)  computation  of  spectral  ratios 
from  long-period  body  wave  components. 

(2)  A  data  archive  of  70  mm  film  chips  from  the  WWSSN 
for  selected  large  earthquakes  and  blasts. 

(3)  A  data  archive  of  digitized  long-period  P  waves 
(3-component)  from  selected  stations  of  the  WWSSN. 

(4)  A  data  archive  of  digitized  long-period  diffrac¬ 
ted  core  waves. 

The  project  objective  to  acquire  and  analyze  a  data 
library  of  digital  long-period  body  wave  signals  from  the 
LRSM  and  LASA  networks  was  never  realized,  owing  to  the 
Impracticability  of  large-scale  access  to  these  collections. 

Graduate  student  research: 

L.  M.  Cathles  was  supported  by  the  contract  for  one 
year.  He  collaborated  as  a  co-lnvestlgator  on  the  detection 
of  the  core  shadow  (Phlnney  &  Cathles,  1969). 
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Research  by  visiting  scientists; 

The  following  visiting  researchers  were  supported  by  the  contract,  and 

contributed  to  the  project  output. 

• 

Shelton  S.  Alexander:  Six  months  as  visiting  fellow  in  1965,  while  on  leave 
from  the  Air  Force.  Contributed  to  papers  on  diffraction  by  the  core¬ 
mantle  boundary. 

Christopher  Chapman:  Six  months  as  visiting  student  in  1968-69.  Carried  out 
research  on  body  wave  calculations,  leading  to  a  Ph.  D.  thesis  at 
Cambridge  in  September,  1969. 

Paul  Richards:  One  month  as  visiting  student  in  1969.  Collaborated  with  Chapman 
and  Phinney  on  the  coupling  of  P  and  SV  waves  in  heterogeneous  media, 
reported  in  California  Institute  of  Technology  Ph.  D.  thesis,  1970. 
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PUBLISHED  PAPERS; 


1964;  Robert  A.  Phinney,  Struature  of  the  Earth's  Crust  from  Spectral 
Behavior  of  Long-Period  Body  Waves ^  J«  Geophys.  Res.,  69;  (14), 
2997-3017. 


ABSTRACT 


Long-period  P  waves  from  distant  earthquakes  have  been  analyzed  from 
seismo^^rams  recorded  at  Albuquerque  and  Bermuda  in  light  of  Haskell's  theory  of 
the  spectral  response  of  a  layered  crust.  By  using  the  ratio  of  the  vertical 
spectrum  to  the  horizontal  component  spectrum,  we  obtain  a  function  which  depends 
on  structure  beneath  the  station.  Because  of  the  poorly  understood  nature  of  the 
signal  which  follows  the  first  P  wave  motion,  the  methods  of  power  spectrum 
analysis  are  applied  and  a  lag  window  selected  to  discriminate  against  long  time 
correlations  within  the  signal.  Corrections  for  the  differing  responses  of  the 
three  components  are  made  by  using. the  power  spectral  matrix  of  calibration  sig¬ 
nals.  A  range  of  crustal  models  has  been  found  which  agrees  with  the  data.  We 
are  restricted  primarily  in  choosing  models  for  the  lower  half  of  the  crust; 
variations  of  structure  in  the  upper  half  do  not  noticeably  affect  the  theoretical 
curves  in  the  period  range  considered  (0.02  to  0.20  cps).  At  Albuquerque  the 
crust  is  about  40  km  thick  and  the  lower  crust  has  velocities  in  the  range  6.6  to 
7.0  km/sec.  The  Mohorovicic  (M)  discontinuity  under  Bermuda  is  12  km  below  sea 
level,  and  the  structure  appears  to  be  a  normal  oceanic  crust  depressed  elasti¬ 
cally  by  the  weight  of  the  volcanics  which  compose  the  island. 
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PUBLISHED  PAPERS  (cont): 

1965:  Robert  A.  Phinney,  Theoretical  Caloulationa  of  the  Spectrum  of  First 

Arrivals  in  Layered  Elastic  Medims,  J.  Geophys.  Res.,  70:  (20)  5107-23. 

ABSTRACT 

A  practical  method  for  the  calcjlation  of  the  spectral  parameters  of 

first-arriving  signals  in  seismology  has  been  the  object  of  much  theoretical  work 

in  recent  years.  The  difficulty  has  been  that  the  first  arrival  usually  behaves 

as  an  imperfectly  trapped  mode.  Mathematically,  it  arises  from  the  contributions 

of  branch  line  integrals  and  complex  poles.  Attempts  to  transform  the  solution 

into  a  generalization  of  the  normal  modes  have  been  a  mathematical  success  only. 

Because  of  the  complexity  of  this  solution,  a  different,  less  elegant  approach  is 

demanded.  A  practical  technique  is  proposed.  By  a  change  of  variable,  the  twice 

-iuR 

transformed  solution  is  separated  into  a  product  of  the  form  f(u)e  .  This 
can  be  integrated  with  respect  to  the  phase  variable  u  ,  using  standard  quad¬ 
rature  methods,  the  real  part  of  u  changing  most  rapidly  along  the  integration 
path.  By  making  the  frequency  complex,  it  is  possible  to  displace  any  singulari¬ 
ties  away  from  the  vicinity  of  the  contour.  This  gives  the  spectrum  of  the  signal 
as  viewed  through  an  exponentially  decaying  time  window,  naking  it  possible  to 
work  with  the  first  arrival  by  itself. 
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PUBLISHED  PAPERS  (cent): 

1966:  Robert  A.  Phinney,  DiscuBBlon  In  A  Di»ou§§ion  on  Rooont  Adoviooo  in  th» 

Toohniqu*  of  Scimia  Recording  and  Analysis,  Proc .  Ro** .  Soc ,  A . ,  2909. 


ABSTRACT 

Hany  of  the  applications  of  body  waves  to  detection  and  identlf icstion 
of  explosions  have  depended  on  the  understanding  of  the  behavior  of  body  wave 
amplitudes  in  a  radially  heterogeneous  earth.  It  is  possible,  by  means  cC  theory 
now  available,  to  compute  wave-theoretical  body  wave  spectral  amplitudes  which 
take  account  of  focusing  and  diffraction  phenomena  not  described  by  ordinary  ‘.'ey 
theory.  The  existence  of  these  phenomena  demands  that  amplitudes  of  body  waves 
be  determined  only  by  careful  study  of  the  spectral  properties  of  the  signals 
and  after  corrections  have  been  made  for  knotm  distortions. 


lUiuiHD  rmiis  (coot) I 


riool  Koport 
AT 

If  StpiMbtr  1970 


IfMi  SmIioo  t.  AIoModor  ond  Aobort  A.  Phinno/t  A  Study  of  tho  Coro^^tantlo 
Somdary  Voimg  f  Uacooo  Oifftoaud  ty  tho  Sca^th'o  Coro,  J.  Gooohyt  . 
too.,  71i  (}«•), 


AiCTAACT 

Mtloo  of  tho  spoctruoo  of  diffroetod  P  movoo  cor.  bo  uood  to  iooloto 
offocts  of  propofotion  along  tho  coro-oontlo  boundary,  and  tho  rato  of  attonua- 
tloo  with  distanco  aa  a  function  of  froquoncy  is  a  uaoful  quantltotlvo  aoasuro 
of  eoro'oantlo  proportlos.  Spoctruos  of  tho  attonuation  coofflclont  for  a  nunbor 
of  patchoa  on  tho  coro'oantlo  boundary  woro  ooaaurod  in  tho  froquoncy  rango  0.01 
to  0.12  cpa,  ualng  data  froo  largo  oarthquakoa  in  difforont  locations.  Although 
tho  gonoral  trond  of  tho  obsorvations  follows  tho  thoorotically  oxpoctod 
froquoncy  dopondonco,  significant  lateral  variations  occur  both  in  lovol  and 
ahapo  of  tho  attonuation  spoctrua,  inplying  lateral  variations  in  coro^oantlo 
proportlos.  In  particular,  resonance  absM^tion  peaks  are  observed  st  certain 
froquoncios  which  osy  bo  a  consoquonco  of  layering  in  tho  deepest  part  of  tho 
SMOtlo.  Tho  results  are  corrolatod  with  geographic  location  in  an  effort  to 
coswonco  to  Mp  tho  lateral  variations  in  coroHaantlo  properties.  In  addition, 
tho  observed  phase  spoctrusw  indicate  tiiat  tho  diffracted  P  waves  are  dispersed 
little  or  none  at  all  for  core-sMntle  boundary  propagation. 
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PUBLISHED  PAPERS  (cont): 

1966:  Robert  A.  Phinney  and  Shelton  S.  Alexander,  P  Mzue  Diffraotion  Thtory 

and  tfce  Struotur*  of  th*  Coro-Hantla  Boundary,  J.  Geophys.  Rea..  71; 
(24),  5959-75. 


ABSTRACT 

The  theory  for  diffraction  of  elastic  waves  in  a  spherical  geometry  is 
reviewed  and  modified  to  apply  to  a  multilayer  geometry.  Wa  have  made  lumerical 
calculations  of  the  attenuation  function  for  P  waves  in  the  shadow  of  the  core 
for  certain  basic  geometries  and  for  models  which  have  a  single  laye  -  at  the 
core-mantle  boundary.  Resonance  phenomena  due  to  the  presence  of  the  layer  are 
discussed,  and  an  attempt  is  made  to  explain  observed  peaks  in  the  attenuation 
spectrum  for  various  regions  of  the  core  boundary.  A  region  between  and  160 
km  thick,  with  a  reduced  shear  velocity  and  an  increased  density,  appears  to  be 
the  most  likely  explanation  of  the  data. 
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PUBLISHED  PAPERS(cont); 

1969;  Robert  A.  Phinney  and  Lawrence  M.  Cathles,  Diffraction  cf  P  by  the 
Core:  A  Study  of  Long-period  Amplitudes  near  the  Edge  of  the  Shadow ^ 

J.  Geophys.  Res.,  74;  (6),  1556-74. 

ABSTRACT 

Theoretical  expressions  for  the  amplitude  of  P  waves  in  the  neighbor¬ 
hood  of  the  core  shadow  boundary  have  been  formulated  and  evaluated  by  computer. 
Comparison  with  2.5-sec  data  of  Sacks  confirms  his  determination  of  the  shadow 
position  at  96.5®.  A  square-root  frequency  scaling  law  is  found  for  the  high- 
frequency  diffraction  pattern  and  permits  the  use  of  the  Fresnel  integral  function 
at  short  periods.  Lcr.g-period  amplitude  data  for  three  paths  have  been  sp^ctrun- 
analyzed  and  compared  against  the  theoretical  amplitude  dropoff  going  into  the 
shadow.  The  agreement  is  about  as  good  as  the  data  will  permit  and  no  indication 
is  found  at  these  periods  of  a  low-amplitude  defocused  region  at  90°  to  96°, 
after  diffraction  has  been  accounted  for. 
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PUBLISHED  PAPERS  (cont): 

1969:  Robert  A.  Phinney  and  Shelton  S.  Alexander,  The  Effect  of  a  Velocity 

Gradient  at  the  Base  of  the  Mantle  on  Diffracted  P  Vaves  in  the 
Shadow^  J.  Geophys.  Res.-,  74:  (20),  4967-71. 

ABSTRACT 

Asymptotic  expressions  for  the  P-wave  solution  near  a  turning  point 
may  be  combined  with  the  boundary  conditions  prevailing  at  the  mantle-core 
boundary  to  obtain  an  asymptotic  expression  for  the  reflection  coefficient  at 
this  boundary  that  preserves  the  diffraction  effects  of  curvature  and  velocity 
gradient.  Increasing  the  velocity  gradient  is  shovm  to  increase  the  decay  and 
coefficient  for  the  diffracted  wave  in  the  shadow.  Consideration  of  recent 
published  data  and  some  new  short-period  data  suggests  that  regional  differences 
in  the  decay  of  this  signal  may  be  due  to  difference  in  the  gradient  at  the  base 
of  the  mantle. 
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PUBLISHED  PAPERS  (cont)! 

1970:  Robert  A.  Phinney,  Reflection  of  Acoustic  Waves  from  a  Continuously 

Varying  Interfacial  Region,  Rev.  Geophys..  8,  517-532,  1970. 


ABSTRACT 

Reflection-transmission  problems  for  plane  waves  can  be  analytically 
solved  in  a  restricted  number  of  special  cases  for  which  asymptotic  connection 
formulas  exist.  The  technique  is  demonstrated  for  the  case  of  a  parabolic 
barrier.  A  method  due  to  Epstein,  which  nakes  use  of  the  properties  of  hyper¬ 
geometric  functions,  may  be  used  to  obtain  solutions  for  a  variety  of  transi¬ 
tional  structures.  This  method  is  taken  over  from  the  field  of  ionospheric 
radio  wave  propagation  and  adapted  to  the  acoustic  problem  without  neglect  of 
the  density  gradient.  The  solution  for  a  monotonic  transitional  zone  is  obtained, 
and  the  behavior  of  the  subcritical  reflection  discussed,  with  reference  to  the 
inner  core  reflection  as  an  example. 
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MANUSCRIPTS  AND  UNPUBLISHED  REPORTS: 


1969:  Christopher  H.  Chapman,  Seismic  Vave  Diffraction  Theory ^  A  dissertation 

for  the  degree  of  Doctor  of  Philosophy,  Christ's  College,  Cambridge. 

ABSTRACT 

The  problem  of  body  wave  propagation  is  considered  for  a  radially  hetero¬ 
geneous  earth.  The  approach  is  based  on  an  expansion  of  the  formal  solution  into 
integrals  representing  "generalized  rays".  The  various  ray  integrals  for  a  homo¬ 
geneous  mantle/homogeneous  cere  model  have  been  studied  for  both  high  frequencies 
(asymptotic  ray  theory)  and  exactly,  by  nsrerical  integration.  The  numerical 
integration  method  involves  complex  integration  in  the  neighborhood  of  the  saddle 
point  representing  the  geometric  ray.  For  the  above  model,  the  signals  P,  PeP, 
diffracted  P,  ?!<?;  diffracted  PKP,  S,  ScS,  SKS,  SnKS,  were  treated. 

Extension  of  these  methods  was  made  to  the  case  of  a  general  earth  model. 
The  various  analytic  approximations  leading  to  ordinary  ray  theory  and  their  ex¬ 
tensions  were  discussed.  Numerical  contour  integration  of  the  ray  integrals  was 
done  as  above.  The  computation  of  the  response  function  for  the  general  earth 
models  is  done  by  a  combination  of  Runge-Kutta  integration  of  the  equations  of 
motion  and  asymptotic  expressions,  for  regions  where  these  are  valid.  The  follow¬ 
ing  phases  were  computed:  P,  PeP,  diffracted  P,  S,  ScS,  diffracted  S,  PKP, 
PKIKP. 


This  thesis  is  being  divided  into  two  or  three  distinct  papers,  which 
are  being  submitted  for  publication.  The  first  has  already  been  submitted  to 
the  Geophysical  Journal. 
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PAPERS  AT  MEETINGS:  (Not  published  elsewhere). 

1968:  Robert  A.  Phinney,  Princeton  University,  Princeton,  New  Jersey. 

Diffraction  Theory  for  PKP.  Annual  meeting  of  the  American 

Geophysical  Union. 

ABSTRACT 

I  have  computed  numerical  results  for  the  wave-theoretical  PKP  signal 
as  a  function  of  frequency  for  an  elastic/fluid  model  of  the  mantle-core  inter¬ 
face.  This  shows  the  kinematic  behavior  of  the  pulse  in  the  neighborhood  of  the 
caustic  and  puts  limits  on  the  frequencies  for  which  ray  theory  is  adequate.  The 
predicted  ampJ itude  decay  into  the  shadow  region  takes  place  over  several  degrees, 
even  at  periods  of  1-2  sec.  Some  preliminary  analyses  of  Icrg-period  data  are 
available  for  the  vicinity  of  the  PKP  caustic  at  141®. 
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PAPERS  AT  MEETINGS;  (not  published  elsewhere) 


1969;  Robert  A.  Phinney,  r.'’inceton,  New  Jersey  and  Christopher  Chapman, 

Cambridge,  England,  Kinenatioa  of  Shadow  Formation.  Annual  meeting 
of  the  American  Geophysical  Union. 


ABSTRACT 

The  increasing  resolution  and  precision  of  modern  seismic  systems  has 
made  it  necessary  to  more  fully  understand  the  behavior  of  body  wave  amplitudes 
in  the  neighborhood  of  cusps  and  endpoints  of  the  travel-time  curve.  These 
arise  when  the  turning  point  of  a  ray  is  sufficiently  close  to  a  region  where 
the  density  and  velocities  do  not  satisfy  the  smoothness  criteria,  of  JWKB  (ray) 
theory.  Most  situations  of  interest  may  be  modeled  by  regions  where  the  wave 
solutions  may  be  shown  to  be  connected  with  the  interface  waves  traveling  along 
the  structures  of  interest;  Analysis  of  a  wide  variety  of  first  and  second  order 
discontinuities  is  possible.  An  example  is  given  showing  the  role  of  the  velo¬ 
city  gradient  in  determining  the  decay  of  the  diffracted  P  phase  in  the  core 
shadow . 
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RELATED  WORK: 

1968:  Michael  Yates,  Washington,  D.  C.  and  Robert  Phinney,  Princeton, 

New  Jersey.  Elastic  Anisotropy  in  Rocks  from  the  Stillwater  Igneous 
Complex^  Montana^  and  the  Tinaquillo  Peridotitej  Venezuela.  Annual 
meeting  of  the  American  Geophysical  Union. 

ABSTRACT 

Longitudinal  seismic  velocities  have  been  measured  in  24  cores  from  8 
sites  in  the  Stillwater  igneous  complex,  Montana,  and  in  3  cores  from  a  sample 
of  the  Tinaquillo  peridotite,  Venezuela.  The  velocity  anisotropy,  AV/V  ,  for 
the  Stillwater  sites  ranges  from  2%  to  5%.  Correlation  with  petrofabric  diagrams 
shows  that  bronzitites  have  minimum  velocities  normal  to  the  layering  in  the  Still¬ 
water  due  to  preferred  orientation  of  the  b  axes.  Anorthosites  show  a  velocity 
maximum  normal  to  the  layering  with  c  crystallographic  axes  concentrated  in  the 
layering  plane.  The  major  characteristics  of  the  preferred  orientation  of  both 
bronzite  and  plagioclase  are  consistent  over  distances  of  14  miles  and  support  a 
gravitational  settling  mechanism  for  the  layering.  Although  the  Tinaquillo  peri¬ 
dotite  samples  appear  to  have  experienced  intensive  shearing  in  their  history,  a 
lack  of  preferred  orientation  in  the  olivine  is  inferred  from  the  low  seismic 
velocity  anisotropy  (2%).  This  is  interpreted  as  being  due  to  a  degradation  of 
shear-induced  preferred  orientations  by  continued  shearing. 
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ABSTRACT 

The  problem  of  determining  the  refractivity  profile  of  a  planetary 
atmosphere  from  optical  or  radio  occultation  data  is  identical  in  principle  to 
the  problem  of  determining  the  variation  of  seismic  velocities  in  the  earth  from 
the  observed  travel  times  of  seismic  body  waves.  In  either  case,  a  complete  set 
of  data  can  be  inverted  uniquely,  the  only  constraints  being  those  fundamental 
to  geometric  optics.  Expressions  are  given  for  converting  observed  Doppler  shifts 
to  the  index  of  refraction  as  a  function  of  depth  in  the  atmosphere.  The  effect 
of  various  approximations  on  the  analysis  is  discussed;  it  is  found  that  a  "thin 
atmosphere"  approximation  simplifies  the  matherratics  and  preserves  the  singularity 
at  the  critical  ray  curvature. 
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13.  ABSTMACT 


This  final  report  describes  research  completed  under  Contract  AF49(638)-1243 ,  during 
the  period  1  July  1903  to  31  December  1969.  The  primary  objectives  of  the  project 
were  as  follows: 

a.  The  development  of  theoretical  methods  for  treating  the  amplitude  and 
phase  characteristics  of  body  waves  in  the  earth, 

b.  The  development  of  parallel  numerical  methods  for  computing  the  amplitudes 
and  pha.sr.s  of  body  waves  in  a  radially  hecerog.jneous  earth. 

c.  Tl.e  applicatloii  of  these  rebulLs  to  a  iiiore  complete  unuetbLandlng  of  obseivcd 
lo.ig  and  short  period  signals  frotn  explosions  and  earthquakes. 


Key  words: 

Long  Period 
Seismogram 
P  waves 

Radially  heterogeneous  earth. 


DD  ,'r..i4-/3 


UNO  AS  'j)  rCD 


